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1. Introduction 

Production cross-sections and decay widths of the Standard Electroweak Model Higgs boson 
are nowadays among the most extensively analyzed theoretical quantities (for recent reviews, see, 
e.g., [[T]], [^). Indeed, the main hope of scientific community is that this essential ingredient of 
the Standard Model may be discovered, if not at Fermilab Tevatron, then at the forthcoming LHC 
experiments at CERN. There is great interest in the "low-mass" region 1 14.5 GeV< Mh < 2Mw, 
because a "low-mass" Higgs boson is heavily favored by Standard Model analysis of the available 
precision data. The lower bound, 114.5 GeV, was obtained from the direct searches of Higgs boson 
at the LEP2 -collider primarily through Higgs boson decay into a bb pair. The decay mode 
H — > bb, which is dominant for the "low-mass" region, is important for Higgs boson searches in 
certain associated (semi-inclusive) Higgs boson production processes at the Tevatron and the LHC. 
It is also the main decay mode for diffractively produced Higgs boson searches in CMS-TOTEM 
and, possibly, FP420 experiments. 

It should be stressed, that the uncertainties in r{H —^ bb), analytically calculated in QCD 
using the MS-scheme at the a^^-level [^, dominate the theoretical uncertainty for the branching 
ratio of H 77 decay, which is considered to be the most important process in searches for a "low 
mass" Higgs boson by CMS and ATLAS collaborations at the LHC. Moreover, since at present 
the QCD corrections to the Higgs boson production cross-sections at Tevatron and LHC are known 
at the next-to-next-to-leading order (NNLO) and, partly, even beyond (see, e.g., it becomes 
important to understand how to estimate the theoretical error-bars of QCD-predictions for both 
the production cross-sections and Higgs-boson decay widths, which are also calculated in QCD 
beyond the next-to-leading order (NLO) level. Here we will focus on the analysis of the concrete 
uncertainties of the QCD predictions for T^^^ = r{H —>■ bb), including those which come from the 
on-shell mass parameterizations of this quantity (previous related discussions see in [||]-[|13]]) and 
from the resummations of the terms, typical of the Minkowskian region (see [|T4[]- [p^]). 

2. QCD expressions for F^j^j^ in the MS-scheme 

2.1 Basic results in terms of running Z7-quark mass 

The QCD prediction for Fjji^^ in the MS-scheme is of the form 



^ Hbb - A -T— 

™b 



l + XAriai(MH) 



(>i 



(2.1) 



Here Fq = 3\/2 /SnGYMum^, Mb and Mh are the pole Z?-quark and Higgs boson masses, ^^ (Mh) = 
0Cs(MH)/?r and inb(MH) are the QCD running parameters, defined in the MS-scheme. The coef- 
ficients AFi can be expressed through the sum of the following contributions: the positive contri- 
butions df , calculated directly in the Euclidean region, and the ones proportional to n^- factors, 
which are typical for the Minkowski time-like region. 

The corresponding expressions were derived at the -level in Ref. [|l^] and have the following 
form: 
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Ar4 



df- 



_ 17 

7o(i8o + 27o)7rV3; 



df - [dE(j8o + Yo){po + 270) + i8i7o + 271 (iSo + 270)] 7i^/3 ; 
d4 - [df (iSo + 7o)(3i8o + 270) + df/3i (5/3o + 670 )/2 
4df 7i (j8o + 7o) + p27o + 271 (jSi + 7i ) + 72 (3j8o + 470)] TrVs 
7o(j8o + 7o) (j8o + 270) m + 270) ^V30 , 



(2.2) 

(2.3) 
(2.4) 



(2.5) 



where the nf -dependence of df (2 < / < 4) was evaluated in [20], [[T^] and (the detailed results 



are presented in [p^l]). The coefficients jS, and 7 are the perturbative coefficients of the QCD 
renormalization group (RG) jS -function and mass anomalous dimension function Ym of the MS -like 
schemes. The QCD j8-function will be considered at the 5-loop level: 



da. 



(2.6) 



The expressions for jSq and jSi are well-known and are scheme-independent. The coefficients ^2 



and J83 were analytically evaluated in p2[ ] and p3[ ] and confirmed by independent calculations at 



the level of 3-loops [g^l and 4-loops [g^. The 5-loop coefficient ^4 is still unknown, and will 
be estimated using the Pade approximation procedure, developed in [|2^]. The mass anomalous 
dimension function is defined as 



dlnfx^ 



7m 



-Jo as 



7i a] ■ 



72 a,^- 73^^-74 «.^ + 0(af) • 



(2.7) 



The 4-loop correction 73(2^ was independently calculated in [ p7| ] and [28]. At 5-loops 74 may 



be modelled using the Pade approximation procedure of Ref.[26] mentioned above. Note that the 



consideration of the explicit nf dependence of the 5-loop coefficients of Eq.(2.6) and Eq.(2.7) pre 



sented in Ref.[13] strongly suggests that the available Pade estimate of the 74-coefficient contains 



more uncertainties than the Pade estimate of the coefficient j84, given in Ref.[g^]. This conclusion 
is supported in part by the fact that the Pade estimated value of the n^-part of 74 [26] is over 3 times 
smaller, than the result of the explicit analytical calculation of Ref. [E^. It should be stressed. 



however, that the uncertainties of the estimated 5-loop contributions to the QCD j3 -function and 
mass anomalous dimension function 7,,, are not so important in the definition of the running of the 
^-quark mass from the pole mass mb to the pole mass of Higgs boson Mh. This effect of running 
is described by the solution of the following RG equation: 



mb(MH) =mb(mb)exp 



dx 



-mg(mb) 



fl'.v(mb 



AD(a, (mb)) J 



(2.8) 



where AD{as) is a polynomial of 4-th order in the QCD expansion parameter = ttv/ K (see [ 13 ]). 



We will complete this section by presenting the numerical values of Eq.(2.2)-Eq.(2.5) in the case of 
nf = 5 active quark flavours and comment on the significance of other QED and QCD contributions 
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to Fjjbij, first considered in Refs.[|8[],[pO|], [^], [^]. In tlie Higgs boson masses region of interest, 
the expression for Eq.(^]l]) may be expressed as 



Hbb 



l + £Ariai(MH) 



i>l 



(2.9) 



V2 
~8k 



GFMHm2(MH) 



1 + 5.667a,(MH) + 29.15a,(MH)^ +41.76a,(MH)^ - 825.7a,(MH) 



Substituting the value ^^(Mh) ~ 0.0366 (which corresponds to a,(MH = 120 GeV) 0.115) into 



Eq.(2.9), and decomposing the coefficients in the Minkowskian series into Euclidean contributions 



and Minkowskian-type ;r^-effects, one can get from the work of Ref . the following numbers 



Hbb 



.m^(MH) 
i =i — 



1 + 0.207 + 0.039 + 0.0020 - 0.0015 



(2.10) 



1 + 0.207 + (0.056 - 0.017) + (0.017 - 0.015) + (0.0063 - 0.0078) 



where the negative numbers in the round brackets come from the effects of analytical continuation. 



Having a careful look at Eq. ( [2.10| ) we may conclude that in the Euclidean region the perturba- 
tive series is well-behaved and the Ti^-contributions typical of the Minkowskian region are also 
decreasing from order to order. However, in view of the strong interplay between these two effects 
in the third and fourth terms, the latter ones are becoming numerically equivalent. This feature 
spoils the convergence of the perturbation series in the Euclidean region. Therefore, to improve 
the precision of the perturbative prediction in the Minkowskian region it seems natural to sum up 



these 71 - terms using the ideas, developed in the 80s (see, e.g., [33], [M], [35] and [36]). Due 



to the works discussed, e.g., in []37[], these ideas now have a more solid theoretical background. 
We will describe some applications of these resummation procedures to r{H bb) later on. Here 
we stress that the truncated perturbative expansions of Eq. (5^91) have some additional uncertainties. 



These include Mh and ^-quark mass dependent QCD [31], [32] and QED [3C] contributions: 



Hbb 



V2 
"87^ 



GFMHm2(MH) 



(2.11) 



where At and A^^^ is defined following Refs. [02|], [|^] as 



At = a; (3.111 -0.667L,) + ^(-10 + 4L( + -/?i(mg/M^)) 



M^ 



(2.12) 



+ a'^ 50.474 - 8.167Lf - 1.278Lf ]+ai—f[ 0.241 - 0.070L, 



mr 



+ XAl -4.913a, + a?(-72.117-20.945L, 



Lt = /?i(M^/m^), Xf = G-prn^ / {Sk^ V2) , mt is the f-quark pole mass, mb = mb(MH) 



m? 



AQ™= 0.472-3.336^ a - 1.455a^ + 1.301 



M2 



(2.13) 
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Using a = a(MH)/7r=0.0027 ( a(MH)"^ w 129), mt = 175 GeV, Mh = 120 GeV, mt = 2.8 GeV, 
Of = 1.1667 X 10"^ GeV"^ 



At 



4.84 X 10 



-3 



1.7 X 10 



-5 



+ 2.27 X 10^^ + 1.85 X 10" 



+ 3.2 X 10 



-3 



5.75 X 10 



-4 



2.42 X 10 



-4 



^QED 



1.1 X 10"^-4.5x 10"^-9x 19"'^-1.2x lO"'* 



(2.14) 

(2.15) 
(2.16) 

(2.17) 



Comparing the numbers presented in Eq.(2.10) and Eq.(2.15)-Eq.(2.17), we conclude that it seems 



more natural to take into account order -terms in Eq.(2.10) only after the possible discovery 



of the Standard Model Higgs boson . Indeed, one can see, that even for the light Higgs boson the 



numerical values of the order -contributions to Eq.( 2.10 ) are comparable with the leading Mh 



and mt- dependent terms in Eqs. ( |2. 14| )-( |2. 16| ) and with the leading QED correction in Eq.( |2. 17| ). 
These terms can be neglected at the current level of the experimental precision of "Higgs-hunting". 

2.2 The relations between different definitions of b-quark mass. 

In the discussions above we used two definitions of the Z^-quark mass, namely the pole mass 
mb and the running mass iiTb. At the maximal order we are interested in, i.e. at the next-to-next-to- 
next-to-leading order (N^LO), these two definitions are related in the following way 

mb(mb) 



1 



^a.v(mb) 



18.556a,(mb)^ - 175.76a,(mb)^ - 1892a,(mb)' 



(2.18) 



where the first three coefficients come from the calculations of Refs. [38], [p^], while the numerical 
estimate of the a^'^-one is the updated variant of the estimate of Ref. [19]. It takes into account the 



explicit expression for the 0{a^) term in Eq.(2. 18) in the effective charges procedure applied in 
Ref.[[l9|] and developed developed previously in Ref. [40]. The pronounced feature of Eq.( 2.18 ) is 



the rapid growth of the coefficients in this relation. This property agrees with the expectations for 
the fast increase of the coefficients of this perturbation series, revealed in the process of applications 



of the QCD renormalon approach in Refs. [41], mM (for a discussion see [43]) 



There are different points of view concerning the application of various definitions of Z>-quark 
mass in phenomenological studies. 

1. The most popular one is that in view of the factorial growth of the coefficients evident in 
Eq.(2.18) it is better to avoid application of the pole mass mb and to use instead the MS- 
scheme running Z?-quark mass mb(;U) normahzed at the scale = mb (see e.g. [44], [p5t]). 



2. 



It is also possible to consider the invariant Z^-quark mass, which is related to the running 
mass, normalized at the scale = mb as 

-1 



mb = mb(mb) 



a,(mb)''oAD(a,(mb)) 



(2.19) 



where AD is defined in Eq.( |2.8| ). The concept of the invariant mass is rather useful in treating 
TT^-contributions to r{H bb) (see [|14|]- [jl|]). 



5 



Uncertainties of QCD predictions for Higgs boson decay into bb 



A .L .Kataev 



3. The pole b-qasak mass is frequently used in the MOM-scheme [|4^]. Within this prescription 
threshold effects of heavy quarks may be understood rather easily [p7|], [p8|], ||49|]. Moreover, 
the concept of the pole Z?-quark mass is commonly applied in considerations of deep-inelastic 
scattering processes [|^] , and what is even more importantly for the LHC, in global fits of 



parton distributions [51 1. Note, however, that quite recently the three-loop transformation 



of the MOM-scheme to the MS-scheme was analysed in Ref.[]52|] where the equivalence of 
these two approaches outside the threshold region was demonstrated. 

4. Keeping in mind the advantages of both running and pole definitions of the b-quaik mass one 
may analyze the effects of the RG resummation of a" In'" -terms {m< n) by comparing 

theoretical predictions for concrete physical quantities, which depend on iiTb and mb. 



Note, that for Eq.( |2.18 ), Eq.( 2.19| ) and other perturbative series, discussed in this work, we coor- 
dinate their truncated expressions with the truncation of the inverse logarithmic expressions for a^ . 
At the N^'^^^^LO of perturbation theory (1 < k < 5) they are defined as 



A)Logi 



iNLO 



A)Log2 



j8iln(Log2) 
i8o'Log2 



)n2lo 



ai(At^)NLO + Aa,(/i^)N2Lo av(jU^)N3LO = ^^(M^ 



)n4lo = <^s{^i )n3lo + Aa^lM )n''LO 



(2.20) 

)N2LO + Aa.(At^)N3LO (2-21) 
(2.22) 



where Logk = ln(/i^/A^), are the values of the QCD scale parameter A^, extracted from 
the experimental data for concrete physical quantities taking into account N^'^^'^LO perturbative 
QCD corrections, which depend on the number of active quark flavours, nf. The definitions of 



Aai(jU )f^(k-i)Lo 3 < k < 5 in Eq.( 2.21 ) and Eg. ( 2.22 ) contain high-order scheme-dependent coeffi 



cients j8i:_i of the QCD jS -function of Eq.(|2.6|) and A^ as well (for details see Ref. []13|]) 



MS 



2.3 Explicit expression for in terms of pole mass 



Consider the MS-scheme perturbative series for T^^^ from Eq.(2.1). To transform it to the 
case when instead of the running mass mb(MH) the pole mass mb is used, one should make the 
following steps: 

• express mb(MH) in terms of mb(mb) and (mb) by solving the RG equation for the running 



mass, defined in Eq.(2.8); 



apply Eq.(2. 18), which relates the square of the running mass mb(mb) to the square of the 
pole mass mb via a perturbative expansion in powers ^^(inb) = o;s(mb)/7r; 

reexpress powers of ^^(mb), which appear at the first and and second steps, in terms of 



powers of fl'.s(MH) using the RG equation of Eq.(|2.6|) for the QCD coupling constant, 



The explicit solutions of equations mentioned above, namely the solutions of Eq.([2.8|) and Eq.(|2.6[), 



were written down in Ref.[]53|] and extended to higher order level in Ref.019|]. Their application 
result in the appearance in the expressions for Tfjbb 



Hbb 



(2.23) 
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of the RG controllable L = In (M^/m^) -terms , which enter into the coefficients AF^ in the follow- 



ing way 



Art 
Ar^ 



3-2L; 

-4.5202 - 18. 139L + 0.08333l2 ; 
-316.88- 133.42L- 
-4366.2- 1094.6L- 55.867 - 1.8065L^ +0.0477L4 . 



1.1551L2 +0.0509 L^: 



(2.24) 
(2.25) 
(2.26) 
(2.27) 



The numerical values of the calculated logarithmic terms are not small. However, in the case of 
i > 2 they tend to cancel each other in the final results for APj'. 

The related variant of Eq.(p]), where the RG-controllable terms are summed up, may be 
written down as 



Hbb 



a,(mb) 



(24/23) 



AD(a, (mb))2 



l + ^Ariai(MH) 



i>l 



(2.28) 



where 



X (l - -a,(mb) - 18.556a,(mb)2 - 175.76a,(mb)^ - 1892a,(mb)^) 



AD{a,f = l+2.351a,+4.383a? + 3.873a?- 15.15a| 



(2.29) 



We will factorize the term Fq out of the expressions for Eq.(2.23) and Eq.(2.28) as well. These 
representations are rather convenient for comparing different parameterizations of ^^^^^^ and of the 
ratio R(MH)=rHbb/rS. 

2.4 Values of the QCD parameters used 



To analyze the behavior of the truncated series in Eq.(2.23) and Eq. (2.28) and of the various 
approximations for the R-ratio it is necessary to fix definite values of the QCD parameters mb, 
A^^"*' and A^"^'. This is done in Table 1. 



order 


mb GeV 


AgL-) MeV 


AgL-) MeV 


LO 


4.74 


220 


168 


NLO 


4.86 


347 


254 


N^LO 


5.02 


331 


242 


N^LO 


5.23 


333 


243 


N^LO 


5.45 


333 


241 



Table 1: The values of the QCD parameters used. 



The LO, N'^LO (1 < k < 3) results for the pole mass mb are taken from Ref. [54], where they 
were extracted using the relation between the mass of the T(1S) - resonance, mb and the ground state 
energy of T(lS)-system. The N^LO estimate of mb is our theoretical guess. The LO, NLO, N^LO 
values for Ag^'*\ given in Table 1, come from the recent parton distribution fits of Ref. [p^]. 



At the NLO and N"^LO level these numbers agree with the values of A^ , which were extracted 

° MS ' 

from the first N^LO QCD analysis of the experimental data, performed in Ref. [56]. The fitted 
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data were the Tevatron experimental data points for the xF^ structure function of the vA'^ deep- 
inelastic scattering (DIS) process and were obtained by CCFR collaboration [57]. The N^LO value 
of A^^"^^ in Table 1 is one of the results of Ref.[l5^. In view of the indications for a convergence 

. (nf=4) 



MS 



of the fits performed at N LO revealed in Ref. [ ]56| ] we will use the same value of Ai^ ' as at 



the N^LO level. To get the results for A^^^\ given in the last column of Table 1, we apply the 
NLO and N^LO matching conditions of Ref. [58] (the latter ones were corrected a bit in Ref. [59]). 
At the N^LO we use the expressions from Ref. [^]. At the N^LO the analytical relation from 



Ref. [61] was applied. Note, that the calculations of Ref. [61] confirmed the validity of analogous 
expressions, obtained (partly numerically) in Ref. [62] and Ref.[^3|]. The related NLO-N'^LO results 



for as(Mz) are contained in the interval (0.118-0.119). Thus they agree with the world average 
value of av(Mz) (for a recent review see Ref.[Q]). 

2.5 The comparisons of the renormalization group improved and the truncated pole-mass 
parameterizations 




Figure 1: The quantities analyzed in the pole (or on-shell) mass approach. 




Figure 2: The quantities analyzed in the approach with explicit RG-resummation 



Different perturbative approximations for two parameterizations of Fj^^^, (see Eq.(2.23) and 
Eq.(2.28)) and of the related ratios R(Mh) are compared in the plots of Figure 1 and Figure 2. 
Looking carefully at these plots, we conclude that 
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1. the partial cancellation of the RG-controUable large contributions to Eqs.(2.25)-(2.27), which 
are proportional to L = ln(M^/m^) and which appear in the coefficients of the pole-mass 
parameterization of F^^i^, leads to a reduction in the size of perturbative corrections. Taking 
them into account results in a decrease of the difference bewteen the behaviour of the curves 
of Fig.l and Fig.2. This feature demonstrates the numerical importance of step-by-step 
application of the RG-resummation approach. 

2. The cancellations between the contributions proportional to L in the coefficient result 
in the small value of the whole ai(MH)'* correction to the perturbative expression for F^^^fj. 
This fact demonstrates the convergence of different theoretical approximants. From a phe- 
nomenological point of view this means that at present the corrections may be neglected 
in the related computer codes for calculating some branching ratios, e.g., for H ^ yy,H ^bb 
and H —>■ T^T^ processes. This conclusion is consistent with our similar statement, made in 
the case of using running Z7-quark nib, but motivated by different theoretical arguments (see 
the end of Sec. 2.1). 

3. The behavior of the RG-resummed expressions for Fj^^i^ and Rfjbb more stable than in the 
case, when RG summation of the mass-dependent terms is not used (see Fig.l). This feature 



supports the application of the complete RG-improved parameterization of Eq.( |2.28| ), which 



is closely related to the one, defined through the running ^-quark mass (see Eq.( |2.9| )). 



4. We observed the existence of a difference AFjj^(j between the truncated pole-mass approach 
and the RG-improved parameterization of F^^j^. A pleasant feature is that this difference 
becomes smaller and smaller in each successive order of perturbation theory considered. 
Indeed, for the phenomenologically interesting value of Higgs boson mass Mh = 120 GeV 
we find that at the a^^-level AFj^,^^ 0.7 MeV, while for the a^^-curves it becomes smaller, 
namely AFj^j^j, ^ 0.3 MeV. At the a^-'-level of the RG-improved MS-scheme series one has 
^Hbb ~ 1-85 MeV for Mh = 120 GeV. For this scale the value of F^^,, with the exphcit 
dependence from the pole-mass is 16 % higher, than its RG-improved estimate. 

We hope, that further studies may clarify whether it is possible to formulate more well-defined 
procedures for determining the theoretical errors for this important characteristic of the Higgs bo- 
son and for the predictions of other perturbative QCD series as well. Indeed, even the MS-scheme 
expression for Fj^,^,^ at the q^^-level reveal additional theoretical uncertainties, which are related not 
only to the the estimates of theoretical errors for the MS-scheme parameters mb(mb) and as{Mz) 
(for a discussion see e.g. Ref.[|65|]). These not previously specified uncertainties result from appli- 
cation of different approaches to the treatment of the typical Minkowskian tt^ -contributions in the 
perturbative expressions for physical quantities. In the case of Fjjgj, the results from application of 
these different procedures will be compared below. 

3. Resummations of ;r^-terms 

3.1 The definitions of the resummed approximants 

As was demonstrated in Ref.[|3|], the kinematic tt^ -terms in the coefficients of perturbative 
series for F^^,^ become comparable with the Euclidean contributions starting from the N^LO 
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corrections (see Eq.(2.1C)). In general, Ti^-terms are resulting from analytical continuation to the 
Minsowskian region of theoretical expressions for physical quantities, defined through the the two- 
point functions of quark or gluon currents, which are calculated in the Euclidean region. These 
terms are starting to manifest themselves from the N^LO. In Ref. |14] the idea of resummation 
of kinematic Ti^-effects of Ref.[35| and Ref. [ |3^ was generalized to the case when the related 
RG-equation has the non-zero anomalous dimension. In that work the case of the N^LO approxi- 
mation for Fjj^i^ was analyzed. However, since the negative -contributions to the N^LO correc- 



tion to rjrTT, turn out to be smaller than the value of corresponding Euclidean term (see Eq.(2.10)), 



the possible development of the tt^ resummation procedure was overlooked by the authors of the 
works of Ref. [14] and Ref. [20], aimed at the study of scheme-dependence of approximations 



for this Higgs boson characteristic. The appearance of Contour Improved Perturbation Theory 



(CIPT) and its application in the semi-hadronic decay channel of the T-lepton []66|], [ ]67| ] pushed 
ahead the real interest in the development of tt^ -resummation approaches both in theoretical and 
phenomenological investigations. 

In order to resum the kinematic tt^ -contributions to perturbative predictions for Fj^i^i^ the au- 
thors of Ref. [|6^ supplemented the CIPT method with the procedure of "Naive Non-Abelianization" 



(NNA) [68], commonly used in the renormalon calculus approach ( a detailed discussion of this 
method can be found in Ref.[|69|]). As the result, the following approximation for ^yg^^^ was obtained 



[15]: 



*,2 r 



-BKM 
Hbb 



2 



(a,(MH))^''Ao + £ (a,(MH))^''dEA„(fl,(MH)) 



n>l 



where the An-functions are defined as: 



[I+P^kW] ^'^'(a,)"-^sin(^arctan(/3o7ra,)) , 



(3.1) 



(3.2) 



Here 5n = n + Vo — 1 and Vo = Ijo/Po depends on the first coefficient of the QCD j3-function 
j3o = (II -2/3nf)/4, introduced inEq.(g||) and ^.(Mh) = l/{l5oln{M^/A^)) It should be stressed 
that the NNA approach is dealing with the leading terms in expansions of perturbative coefficients 
for Euclidean quantities in powers of number of flavor nf and it provides the basis of "large jSq- 
approximation". Within this approximation it is assumed that the terms, proportional to jSq ( where 
p is "large", namely 1 < p < oo) give a qualitatively good approximation for the structure of the 
Euclidean perturbative contributions to physical quantities under study. Indeed, as was shown in 



Ref. [ [70| ] this approach gives correct both for sign for order of magnitude estimates of the per- 
turbative coefficients dp to various physical quantities, including d4- contribution to Fj^^b' defined 



in Eq.(2^). The explicit calculations of Ref.[^] demonstrate, that at nf=5 the real value of 



coefficient is higher, than its NNA estimate from Ref.[|70|] by the factor 5 approximately. ^. 

Within this "large /3o-approximation" is seems more consistent to approximate the perturbative 
QCD expansion parameter ^^(Mh) = o;s(MH)/7r by its LO-expression of Eq.(2.20). Fixing now 
n = in Eq.(3.2) and expanding Aq to the first order in a^, the authors of Ref.[15 ] got 



Ao 



1 



sin(b() arctan(7r/L]viH)) 



boLMH 



bo 



:i + 7rVLMH')''<'/2 



(3.3) 



Unfortunately, NNA results from Table I of the important work of Ref. [H] contain misprints. 
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where bo = Vo — 1, Lmh = ^^i^u/^^)- ^^^^ concrete expression was first derived in Ref. [|T^]. 
Otlier ways of resumming tt^ -contributions to ^^^^^^ were considered within the framework of 



Fractional Analytical Perturbation Theory (FAPT) [ ]16[ ] and its variant (for a review of FAFT see 



Ref. [ |17[ ] where "flavor-coiTccted global FAPT" was also proposed). It should be stressed that the 
cornerstone of FAPT is the Analytical Perturbation Theory approach, which was developed in the 
studies, initiated by the work [f/ll]. 

Applying the "large jSq -expansion", which is equivalent to the choice of the LO expression for 



one can get the FAPT analog of Eq.(3. 1 ). It can be written down as 



p(l;FAPT) _pb'^b 



n>l 



^Svo(Mh) 
71" 



(3.4) 



This approximant almost coincides with the expression of Eq.(3.1). This conclusion was made 



in Ref. [16] taking into account the following definitions of the functions in Eq.(3.4), namely 



2li'^(MH) 
2i2v(Mh) 



sin[(v- l)arccos(LMH/-v/7r2 + L2j^)] 



;r(v-l)(7r2 + L2^jv-i/2 

r(v) f d 



r(?l + v)V dL; 



(3.5) 
(3.6) 



The FAPT N^LO expression of Ref.(3.4) was given in Refs.[[T^, [17]. It has the following 



form: 



p(3;FAPT) 
Hbb 



«,2 



,mi 



mt 



<'(MH) + IdE 

n>l 

l(3) 



^(3) 



(Mh) 



71" 



(3.7) 



Within the framework of FAPT the functions *Bjj^yg(MH) (0 < n < 3) absorb the evolution of 
the MS-scheme running b-quark mass and tt^ -contributions in Eq.(^[T|), which are proportional to 
high-order coefficients of RG-functions j8(av) and Y,„{as) in Eqs.(2.1)-(2.5). 

In the "flavour-corrected (f-c) global FAPT", proposed in Ref.[|l^], the FAPT expression for 
Fjj^i^ is expanded into the functions *b{.'] (s), which absorb all nf dependence from the Euclidean 
coefficients d^ in Eqs.(2.3)-(2.5).The explicit expression for this approximant reads [18] 



-,(3;FAPT,f-c) 
Hbb 



ml 



«i?(mh)+i: 



3 <B 



(3) 

n+Vo;dE 



(Mh) 



ii>l 



;r" 



(3.8) 



The more detailed study of the characteristic feature of formalizm [p^], [ [T8| ] is now in progress 



[72] 



3.2 The comparison of the results of Ti^-resummations. 

We now consider the result of applying these different procedures for resummations of kine- 
matic TT^-effects in the perturbative coefficients for F^i^,^. We will follow the results obtained in 



Ref. [ 18] and discuss the comparison of the behavior of various parameterizations for this quantity, 
which were defined in the previous section, with the aj- and a^- truncated MS-scheme expression 
of Eq.(2.1). 
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In the process of these discussions we use plots, similar to those, presented in Ref.[|T8p, but 
modified at our request by changing the values of the QCD parameters to the ones presented in 



Table 1 . The results of these additional studies were kindly communicated [ l72| ] to us in the form 
of the figures presented below. 







3 

2 . 5 
2 


r„^um) |McV] 

^'-""'^ 
.-^ ^-""''^ 


.'''^ 






1 . 5 


--^ ^^■''■^ 
.-^ ^-"^ 
-■^ ^--"'^ 














Mu [GcV] 


1 




Mil [GcVJ 



IDD 12D 14 



160 180 



80 100 120 140 160 18D 



Figure 3: The comparison of N^LO approximations for "^yvo^- On both panels the dashed Hne is the truncated 
a^^-approximation of the MS-scheme result. The dotted line is the 1-loop FAPT expression of Eq.(3.4), 
which is identical to the CIPT expression of Eq.(3.1). The solid line on the left figure displays the 2-loop 
FAPT analog of Eq.(3.7), while on the right figure it corresponds to 2-loop variant of the "flavour-corrected 
global FAPT" approximant of Eq.(3.8). 



3 


rH^»(Af„) [MeV] 


3 


r„^a(WH) |McV] 


2.5 




2.5 




2 




2 




1.5 




1 . 5 




1 


.1/,., [GoV] 

80 IDD 12D 140 160 180 


1 





Figure 4: The N LO version of Fig. 3. The dashed line is the truncated MS-scheme expression. 

A careful look to the curves from Fig. 3 and Fig.4 leads us to the following observations: 

1. at the N^LO-level the resummation of the Ti^-terms do not lead to detectable effects. The 
related curves almost coincide with the MS-scheme approximation for T^^, obtained in 



Ref. [gO|l; 

2. at the N-^LO level the effects of resummation of the kinematic tt^ -contributions are visible. 
Indeed, for Mh = 120 GeV the value of the resummed approximant is over 0.2 MeV lower. 



than the value of the truncated a?^ MS-scheme expression for Fjj^i,, obtained in Ref. 1 21]. The 
behavior of the curves of the left plot from Fig.3 indicate that for a Higgs boson with this 
mass, the value of the -approximation for its decay width to ^ ^ quarks is ^y^^ « 1 .7 MeV, 
while the resummation of the Ti^-effects decreases this by over 11%; 



3. an obvious message, which comes from the consideration of the right-hand plots of Fig.3 and 
Fig.4 is that the application of "the flavor-corrected global FAPT" with full analytization of 
nf dependence in Tyvoh l^^^s to smaller values of the decay width for Higgs boson, than in the 



case of the truncated MS-scheme approach, and 2-loop and 3-loop FAPT approach [16]. An 
attempt to explain the phenomenological reason for this reduction is now in progress [^]. It 
may be also of interest to compare FAPT applications with the existing applications to Ty^^^ 



of the CIPT based resummation procedure from Ref. [73]. 
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Conclusions. 

In the discussions presented above we found that even for the fixed values of the QCD pa- 
rameters used in Table 1, which do not take into account the existing theoretical and experimental 
uncertainties in the values of av(Mz) (and thus A^^^\ and of the running and pole Z>-quark 
masses, different parameterizations of the decay width of the H ^ bb process deviate from the 
MS -scheme prediction. In our view, this feature demonstrates the existence of additional theoret- 
ical QCD uncertainties, which are not usually considered in phenomenological studies. We still 
do not know what is the role of similar QCD uncertainties in the theoretical predictions for other 
characteristics of Higgs boson production and its other decays (see however, the recent studies of 
resummation of terms proportional to in the N^LO perturbative QCD predictions for the Higgs 
-boson production cross-section and for the hadronic decay rate r(H gg) [74]). It will also be 
interesting to understand whether these uncertainties survive in the branching ratios of the Standard 
Model Higgs boson for different values of its mass. 
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